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ABSTRACT 

In transmission-mode terahertz time-domain spectroscopy (THz-TDS), the thickness of a sample is a critical 
factor that determines an amount of the interaction between terahertz waves and bulk material. If the interaction 
length is too small, a change in the transmitted signal is overwhelmed by fluctuations and noise in the system. In 
this case, the sample can no longer be detected. This article presents a criterion to determine the lower thickness 
boundary of a free-standing film that can still be detectable by free-space transmission-mode THz-TDS. The 
rigorous analysis yields a simple proportional relation between the sample optical length and the system SNR. The 
proposed criterion can help to decide whether an alternative terahertz thin-film sensing modality is necessary. 
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1. INTRODUCTION 

Teraherz time-domain spectroscopy (THz-TDS) is an efficient tool for characterizing materials with broadband 
coherent terahertz radiation.l^ A femtosecond optical laser is utilised in THz-TDS. Interaction of optical pulses 
with an emitter results in a burst of subpicosecond pulses spanning a range from a few hundred gigahertz to a few 
terahertz. The emitted pulses induce a local change in the electric field at the detector. A pump-probe detection 
scheme resolves the amplitude and phase of a terahertz pulse with a high SNR. In a transmission measurement, 
the generated terahertz pulses are changed in the amplitude and phase after interacting with a sample. These 
changes can be used to extract the complex optical constants, or related quantities, of the sample. 

For free-space transmission-mode THz-TDS, each material has an optimal thickness that yields the lowest 
measurement uncertainty in the optical constants.f^l The noise floor clearly deflnes the maximum thickness of 
a sample that still results in the valid optical parameters.'^ For a thin sample, the measurement accuracy 
significantly reduces, because the signal variation due to fluctuations and noise in the system becomes relatively 
strong compared with the change in the amplitude and phase introduced by the sample. The limitation of 
thin-film characterization exists where the film is so thin that the extracted parameters are no longer reliable. 
Likewise, detecting the presence of a thin film reaches a limit when the sample is no longer distinguishable from 
the background material. Despite that, there is no explicit criterion for determining a lower boundary for the 
thickness of a sample under free-space THz-TDS measurement. 

This article presents a mathematical analysis on the limitation of free-space THz-TDS for thin-film detection. 
Initially, the thin-film condition is expressed in terms of optical constants and their standard deviations, which 
are caused by a variation in the signal amplitude. These standard deviations are then further expanded to yield 
the criterion of the detection limitation. It defines the minimal film thickness and refractive index as a function 
of the system SNR. The derivable condition is intuitive and useful as a rule of thumb in estimating the capability 
of a THz-TDS system. It also suggests when other advanced measurement techniques are necessary to detect 
particular samples. 
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2. LIMITATION FOR THIN-FILM DETECTION 



In the context of the uncertainty analysis, the random error, which encompasses fluctuations and noise, is used 
to evaluate the confidence interval in a measurement. This confidence interval localizes or bounds the expected 
value, i.e., the value that is free from the random error. Therefore, in the sample detection scheme, if the 
confidence intervals of the sample and reference measurements are well separated, it can be inferred that the 
sample is unambiguously detected. Otherwise, if the two intervals get overlapped, it is possible that the two 
measurements share an identical expected value. Hence, in this case, the presence of the sample cannot be 
judged. Fig. [T] schematically illustrates both of the situations. It is worth nothing that the systematic error is 
irrelevant in the detection scheme, since it biases the sample and reference measurements equally. 

Mathematically, the confidence interval for any variable x is defined afl 
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fix < X + 



N 



(1) 



where x = X]il=i ^i/^ is the arithmetic mean from N measurements, Sx is the empirical standard deviation, and 
fix is the expected value. From Eq. [1] and the discussion in the previous paragraph, it can be established that 
the sample with optical constants ns{uj) +jKs{uj) can be distinguished from the reference material with hq + jkq 
by terahertz waves if any of these two criteria is satisfied: 
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where An{uj) = ns{ui) - no and Ak{u) = Ks{uj) - kq; '^ln,K},E,^^^ *{n,«;},_E„t' ^^'^ ^}n,K},i ^^'^ empirical 
variances of the sample optical constants caused by the sample and reference amplitude variations, and thickness 
variation, respectively; Ne^^^ and Ne^^i are the numbers of measurements for the sample and reference signals, 
respectively; and Ni is the number of sample thickness measurements. Eq. [2] means that the difference between 
the sample refractive index or extinction coefficient and that of the reference material must be larger than its 
combined standard deviation, which is mainly contributed by two sources of random error, i.e., signal noise and 
thickness measurement.'^ 
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(a) (b) 
Figure 1. Confidence intervals in the detection scheme, (a) The confidence intervals of the sample and reference measure- 
ments are separated. The presence of the sample is confirmed, (b) The two confidence intervals partially overlap. The 
presence of the sample is unclear. 



of sample and reference measurements are equal. Hence, from the given assumptions, Eq. [2]can be simplified to 



AnH > ^4^, (3a) 



N 

A«H > (3b) 



where S{„^^j,(w) + sfn, and N = TVsam = ^rcf- Typically, it can be shown that 

An{uj) > Ak{uj) , and (4a) 
Sn{uj) < Si^{uj) . (4b) 

The conditions in Eq. |4] implies that the criterion in Eq [3a| is less strict than Eq I3bl As discussed earlier, to 
detect the presence of the sample, satisiying only one criterion, here Eq. [3al is enough. 

From the relation between the standard deviation of the refractive index and that of the reference phase in 
Eq. I15al the criterion for sample detection given in Eq. [3a] can be expressed as 

'^An{uj) > y^Sarg(£;,,f)(t^) . (5a) 

Here, I is the sample thickness, and Sarg(£;rof ) ('-^) standard deviation in the phase of the reference spectrum. 

This equation is the ultimate criterion for any free-standing sample to be detectable by normal transmission- 
mode THz-TDS. It means that a phase change introduced to the terahertz pulse by a sample must be larger 
than the scaled standard deviation of the reference phase. This standard deviation Sarg(£;„t) ('*^) t)e obtained 
from a set of repeated time-domain measurements by using Eq. l9bl or the Monte Carlo method. Increasing the 
number of measurements N does not reduce the minimum sample thickness, but on the contrary accentuating 
Sarg(£;,,f)(w) due to the long-term laser instability.'^J 



3. CONCLUSION 

In this article, a limitation in thin-film detection using ordinary transmission-mode THz-TDS is derived on the 
basis of the uncertainty analysis. Provided that the standard deviation of the terahertz signal is known in priori, 
the proposed criterion can be used to estimate the minimal thickness and refractive index of a free-standing 
film that can be detected with a free-space THz-TDS measurement. In future work, a limitation in thin- film 
characterization and a limitation in the case of substrate-based films will be analyzed. The proposed criterion 
will be validated through experiments with different types of films. 



APPENDIX A. STANDARD DEVIATION IN OPTICAL CONSTANTS 

This section shows a derivation of the relation between the variance in the optical constants and the variance 
in the terahertz magnitude and phase. From the amplitude in the time domain, the variance propagates to the 
frequency domain via Fourier transform. The amplitude and phase variances of the reference and sample spectra 
then produce the variance of the transfer function. Eventually, the variance appears at the optical constants. 
The following derivation is adapted from the existing detailed analysi^ 

The discrete Fourier transform of a time-resolved signal, E{k), is given as 

E{u) = ^i;(fc)exp(-jc^fcT) , (6) 

k 

where k is the temporal index, r is the sampling interval, and E{lu) = Et-{lj) + jEi(io). Assuming that the 
amplitude at each time sample is statistically independent from the amplitude at other time samples, the variances 



of the real and imaginary parts s'^ (w) and s|;.(a;), and their covariance s^;^^;. (w) are, respectively)^ 

Se^{uj) = ^ cos^(wA:T)s|(fc) , (7a) 

k 

4.(^) - ^sin2(c.fcr)4(fc) , (7b) 

k 

SErEii^) = — sin(cjfcT) cos(cjfcT)s|;(fc) , (7c) 



where s|;(fc) is the variance of the time-domain signal E{k). The magnitude and phase of the terahertz signal, 
determined from the real and imaginary parts of the complex spectrum, are, respectively. 



arg(i?(w)) — arctan(i?i(w)/_Ei.(i^)) . 
The variances in the magnitude and phase are therefore given by, respectively. 



'In 



'arg(£;) 



\E\(^) = p^^^3*'[^(w)exp(jwfcr)]4(/c) , 
H = ^^^^5'[i?(^)exp(jwfcr)]4(fc), 
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where SR^ and 9^ denote the square of the real and of imaginary parts, respectively. The transfer function of a 
sample is calculated from the reference and sample complex spectra as follows 



ln\H{Lu)\ = ln|i;3am(c^)|-ln|-B,.cf(w)| , 
arg(iJ(w)) = arg(£'sam('^)) - arg(i;rof(cj)) 
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The variances of Eq. llOal and llObI are, respectively, 

sfn\H\{^) = S^n\E,^„,\i^) + Sln\E,,t\i^) ^ (Ha) 

If the reference and sample measurements have a comparable amplitude. Eg. Illal and lllbl can be simplified to 
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The refractive index and the extinction coefficient are evaluated from the transfer function via 
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Thus, the variances of the refractive index and the extinction coefficient due to the magnitude and phase variances 
are 



(^) 4g(H)M, 
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By substituting Eq.[T2]into[T31 the amplitude- related standard deviation in the optical constants can be expressed 
as 



s„(iu) 



'arj 
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Eg 1151 concludes the derivation of the standard deviation in the optical constants. 
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